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Summary. — MCSTHAR++ is a new Monte Carlo code implementing the Statistical 
Hadronization Model. This model assumes that hadronization proceeds through 

the microcanonical decay of massive extended clusters. Unlike other hadronization 
models, in this approach very few free parameters are needed, as has been demon- 
strated in previous studies. The tuning of the model and the comparison with the 
data is ongoing. 

PACS 12.40.Ee - Statistical models of strong interactions. 

PACS 13. 66. Be - Hadrons production by electron-positron collisions. 



1. — Introduction 

The transition from the partonic state to the confined state is described using various 
phenomenological models Hke the Cluster Models [1][2] and the String Model [3]. A 
completely independent model, based on a statistical formulation and known as Statistical 
Hadronization Model [4], has been studied in the literature but it is not available in any 
official release of the Monte Carlo event generators for High Energy Physics. In this paper 
a short overview of the microcanonical formulation of the Statistical Hadronization Model 
will be given and a new Monte Carlo code performing the hadronization process according 
to that model, MCSTHAR++, and its interface to the event generator HERWIG6.510 [5] will 
be described, showing also some preliminary results obtained for the hadronization of 
light quarks only. 
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2. — The Statistical Microcanonical Hadronization Model: a short overview 

As described in [4], the statistical model is based on the hypothesis that in a high 
energy collision some extended objects made of pre-hadronic matter, called clusters or 
fireballs, are produced. These objects are supposed to be colorless and characterized by 
well defined physical quantities like energy, momentum and abelian charges. Each one of 
these clusters hadronizes according to a pure statistical law, which for the microcanonical 
formulation is the following: every multihadronic state confined inside the cluster and 
conserving all the physical quantities of the cluster itself is equally likely. In this picture 
of the hadronization process the cluster is treated as a microcanonical ensemble, whose 
partition function O is given by [6] 



is the sum over all the multihadronic states | hy) confined inside the cluster, while the 
sum on the left equation is over all possible channels, corresponding to the asymptotic 
states, which arc identified by the if-tuple of integer (Nj) = (Ni, N2, Nk), one for 
each hadron species j, and where Pq is the projector on the conserved set of charges. 

The above partition function is then modified introducing the strangeness suppression 
parameter 7^ and multiplying the weight of each channel by the factor 7^= , where Ng 
is the total number of strange and antistrange quarks contained in the hadrons of the 
channel itself. At this point this hadronization model requires only two parameters, 
which need to be tuned on the experimental data: the 7s parameter and the energy 
density of the clusters p, which is taken to be equal for each cluster and which is used to 
convert the mass of a cluster into its volume, a quantity which appears in fl, as can be 
seen in [6]. 

3. — MCSTHAR++: code description and preliminary results 

MCSTHAR++ (Monte Carlo STatistical HAdronization in high energy Reactions) is a 
Monte Carlo code implementing the Statistical Hadronization Model in its microcanoni- 
cal formulation. It is written in Object Oriented C++ and it is built to take as input a set 
of clusters and to give in output a set of hadrons (stable and unstable) , coming from the 
microcanonical hadronization of each one of the incoming clusters. For the hadronization 
of these objects the exact conservation of energy-momentum, electric charge, strangeness, 
baryonic number, charm and beauty charge is imposed. 

In the present case, MCSTHAR++ is interfaced to HERWIG6.510 [5], the clusters used as 
input for MCSTHAR++ are the "standard" HERWIG's clusters (allowed to have non zero bary- 
onic number) and the primary hadrons produced during the hadronization are decayed 
using the HERWIG's routines for the strong and electroweak hadron decay. 

Even if a full tuning of the model is mandatory, it is already possible to see the 
good behavior of the hadronization code in the preliminary results shown in this sec- 
tion. The theoretical predictions are obtained with no tuning at all of the HERWIG's 
parameters and with a reasonable choice of the MCSTHAR++'s parameters: 7s = 0.65 
and p = 0.35 GeV/fm^. The distributions and the multiplicity table show a comparison 
among MCSTHAR++, HERWIG6.510 and LEP (OPAL [7] and DELPHI [8]) data for the 
process e+e~ l/Z^ — ^ uu, dd, ss at 91.2 GeV center of mass energy. 
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Fig. 1. - Charged particle scaled momentum Xp = 2p/^/s and log{l/xp) distribution. 



Table I. - Mean values of the charged particle multiplicity and of the production rate of charged 
pions, charged kaons and (anti)protons, for the hadronization of light quarks only. 





^ch 




Nk± 




MCSTHAR++ 


19.53 ±0.14 


16.64 ±0.11 


1.65 ±0.04 


0.98 ±0.07 


HERWIG6 . 510 


18.601 ±0.006 


15.022 ±0.006 


1.628 ±0.002 


1.736 ±0.002 


DELPHI 


19.94 ±0.34 


16.84 ±0.87 


2.02 ±0.07 


1.07 ±0.05 



4. Conclusions 

MCSTHAR++, a new Monte Carlo code for the Statistical Microcanonical Hadronization, 
has been presented. This code, written in C++, will be soon available to be used with the 
Monte Carlo event generators alternatively to the standard hadronization modules, as it 
has been shown in this paper, where MCSTHAR++ has been interfaced to HERWIG6.510. 

The next steps of this work will be focused on the tuning of MCSTHAR++, completed 
with the hadronization of c and h quarks, on LEP and Tevatron data, and on the com- 
parison with other available hadronization models. 

* * * 
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